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CONSTRAINING THE LIFETIME OF QUASARS FROM THEIR SPATIAL CLUSTERING 
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ABSTRACT 

The hfetime of the luminous phase of quasars is constrained by current observations to be 
10^ ^ tq < 10^ years, but is otherwise unkown. We model the quasar luminosity function in detail in 
the optical and X-ray bands using the Press-Schechter formalism, and show that the expected clustering 
of quasars depends strongly on their assumed lifetime tg. We quantify this dependence, and find that 
existing measurements of the correlation length of quasars are consistent with the range 10^ ^ tq ^ 10*. 
We then show that future measurements of the power spectrum of quasars out to z ~ 3, from the 2dF 
or Sloan Digital Sky Survey, can significantly improve this constraint, and in principle allow a precise 
determination of tg. We estimate the systematic errors introduced by uncertainties in the modeling of 
the quasar-halo relationship, as well as by the possible existence of obscured quasars. 

Subject headings: cosmology: theory ~ cosmology: observation - quasars: formation - large scale 
structure 



1. INTRODUCTION 

A long outstanding problem in cosmology is the syn- 
chronized evolution of the quasar population over the red- 
shift range ^ z ^ 5. Observations in the optical (Pei 
1995) and radio (Shaver et al. 1994) show a pronounced 
peak in the abundance of bright quasars at z ~ 2.5; re- 
cent X-ray observations (Miyaji et al. 2000) confirm the 
rapid rise from z = towards z « 2, but have not shown 
evidence for a decline at still higher rcdshifts. Individ- 
ual quasars are widely understood to consist of supermas- 
sive black holes (BHs) powered by accretion (Lynden-Bell 
1967; Rees 1984). A plausible timescale for quasar activity 
is then the Eddington time, 4x 10^ (e/0.1) yr, the e-folding 
time for the growth of a BH accreting mass at a rate M, 
while shining at the Eddington luminosity with a radia- 
tive efficiency oi L — L^dd — ^Mc^. The lifetime tq of 
the luminous phase of quasars can be estimated directly, 
by considering the space density of quasars and galaxies. 
At z ~ 2, the ratio nq/nc ^ 3 x 10^'^ implies the reassur- 
ingly close value of tq ~ tnuhnq/nq ~ 10^ yr (Blandford 
1999 and references therein). These lifetimes are signifi- 
cantly shorter than the Hubble time, suggesting that the 
quasar population evolves on cosmic time-scales by some 
mechanism other than local accretion physics near the BH. 

It is tempting to identify quasars with halos condensing 
in a cold dark matter (CDM) dominated universe, as the 
halo population naturally evolves on cosmic time-scales 
(Efstathiou & Rees 1988; Haiman & Loeb 1998; Kauff- 
mann & Haehnelt 2000). Furthermore, quasars reside in 
a subset of all galaxies, while the redshift-evolution of the 
galaxy population as a whole (qualitatively similar to that 
of bright quasars) has been successfully described by asso- 



ciating galaxies with dark halos (e.g. Lacey & Cole 1993; 
Kauffmann & White 1993). A further link between galax- 
ies and quasars comes from the recent detection, and mea- 
surements of the masses of massive BHs at the centers of 
nearby galaxies (Magorrian et al. 1998; van der Marel 
1999). 

These arguments suggest that the evolution of the 
quasar population can indeed be described by "semi- 
analytic" models, associating quasars with dark matter 
halos. In this type of modeling, the quasar lifetime plays 
an important role. The quasar phase in a single halo 
could last longer (tq ~ lO^yr), with correspondingly small 
A/bh/-Mhaio ratios, or last shorter (tq ~ lO^yr), with larger 
BH formation efficiencies (Haiman & Loeb 1998; Haehnelt 
et al. 1998). Note that although recent studies have es- 
tablished a correlation between the bulge mass Mbuigc and 
BH mass Afbh, this correlation leaves a considerable uncer- 
tainty in the relation between Mbh and the mass Mhaio of 
its host halo. If the initial density fluctuations are Gaus- 
sian with a CDM power spectrum, the clustering of col- 
lapsed halos is a function of their mass - rarer, more mas- 
sive halos cluster more strongly (Kaiser 1984; Mo & White 
1996). Hence, measurements of quasar clustering are a 
potentially useful probe of both BH formation efficiencies 
and quasar lifetimes (La Franca et al. 1998, Haehnelt et 
al. 1998). 

In this paper, we assess the feasibility of breaking the 
above degeneracy, and inferring quasar lifetimes, from the 
statistics of clustering that will be available from the Sloan 
Digital Sky Survey (SDSS, Gunn & Weinberg 1995) and 
Anglo- Australian Telescope Two-Degree-Field (2dF, Boyle 
et al. 1999). Previous works (e.g. Stephens et al. 1997; 
Sabbey et al. 1999) have yielded estimates suggesting 
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that quasars are clustered more strongly than galaxies. 
However, the current uncertainties are large, especially at 
higher redshifts, where clustering has been found to de- 
crease (lovino & Shaver 1988; lovino et al. 1991), to stay 
constant (Andreani & Cristiani 1992; Croom & Shanks 
1996), or to increase with redshift (La Franca et al. 1998). 
As a result, no strong constraints on the life-time can be 
obtained yet. The key advance of forthcoming surveys over 
previous efforts is two-fold. Because of their sheer size, i.e. 
the large number of quasars covering a large fraction of the 
sky, both shot-noise and sample variance can be beaten 
down, significantly reducing the statistical uncertainties. 
Furthermore, the large sample-size will eliminate the need 
to combine data from different surveys with different se- 
lection criteria, hence allowing cleaner interpretation. 

Recent measurements of the local massive black hole 
density have stimulated discussions of a radiative efficiency 
which is much lower than the usual ^ 0.1 (e.g. Hachnclt 
et al. 1999). A convincing constraint on the lifetime of 
quasars could be therefore highly interesting, as this might 
have implications for the local accretion physics near the 
BH. 

This paper is organized as follows. In § we summa- 
rize our models for the quasar luminosity function, and in 
§ ^ we compute the quasar correlation function in these 
models. In § ^ we compare the model predictions with 
presently available data, and in § [s] we assess the ability 
of future optical redshift surveys to discriminate between 
the various models. In § ^, we repeat our analysis in the 
soft X-ray band, and examine the contribution of quasars 
to the X-ray background, and its auto-correlation. In § ^ 
we address some of the caveats arising from our assump- 
tions, and in § ^ we summarize our conclusions and the 
implications of this work. 

2. MODELS FOR THE QUASAR LUMINOSITY FUNCTION 

In this section, we briefly summarize our model for the 
luminosity function (LF) of quasars, based on associat- 
ing quasar BHs with dark halos. Our treatment is similar 
to previous works (Haiman & Loeb 1998; Haehnelt et al. 
1998), although differs in some of the details. A more ex- 
tensive treatment is provided in the Appendix. The main 
assumption is that there is, on average, a direct monotonic 
relation between halo mass Mhaio and average quasar lu- 
minosity Lm,z, which we parameterize using the simple 
power-law ansatz: 

Lm,z = xo{z)Mhaio (^—^^^ ■ (1) 

Here xo{z) and a{z) are "free functions", whose 
values are found by the requirement that the re- 
sulting luminosity function agrees with observations. 
As explained in the Appendix, our model has one 
free parameter, the lifetime tq, which uniquely deter- 
mines xo{z) and a{z) in any given background cos- 
mology. We assume the background cosmology to 
be either flat (ACDM) with (J^a, f^m, /i, o-g^-i , ") = 
(0.7,0.3,0.65,1.0,1.0) or open (OCDM) with 
(f7A,17ni,/i,cr8h-i,n) = (0,0.3,0.65,0.82,1.3). In 
LCDM, we find (- log[2;o/LoM0^], a) « (1,0.4) and 



w (-0.2,-0.1) for lifetimes of tg = 10* and tg = 
10^'^yr, respectively. Similarly, in OCDM, we find 
(-log[xo/LQM-i],a) « (1,0.25) and « (-0.2,-0.25) 
for these two lifetimes. 

In Figure we demonstrate the agreement between the 
LF computed in our models with the observational data at 
two different redshifts, z — 2 and z = 3. For reference, the 
upper labels in this figure show the apparent magnitudes 
in the SDSS g' band, assuming that the intrinsic quasar 
spectrum is the same as the mean spectrum in the Elvis 
et al. (1994) quasar sample. The photometric detection 
threshold^ of SDSS is g' ~ 22.6, corresponding to a BH 
mass of « 1O*M0 at z = 3 and a three times smaller mass 
at z = 2. As the figure shows, the overall quality of the fits 
is excellent; for reference, the dashed lines show the ad-hoc 
empirical fitting formulae from Pei (1995). Similarly accu- 
rate match to the quasar LF is achieved at different red- 
shifts, and in the models assuming an OCDM cosmology. 
Figure ^ shows, in particular, that the fits obtained from 
the power-law ansatz adopting either a short (solid lines) 
or a long (dotted lines) lifetime are nearly indistinguish- 
able; hence modeling the LF by itself does not constrain 
the quasar lifetime within the limits 10®'^ yr ^ tq <, 10^ 

yr- 

Before considering constraints on the lifetime from clus- 
tering, it is useful to point out that estimates for both 
upper and lower limits on tq follow from the observed lu- 
minosity function alone. 

Lower limit on tq. A halo of mass ^Ihaio is un- 
likely to harbor a BH more massive than w 6 x 
10-3(17fc/r!o)Mhaio = 6 X lO-^Afhaio, whcrc w 6 X 10-^ 
is the ratio Afbh/^^buige found in nearby galaxies (Magor- 
rian et al. 1998) because Mbuige cannot be larger than 
(r2b/r2o)Mhaio- This maximal BH could at best emit 
« 10% of the Eddington luminosity in the B-band, im- 
plying i/Afhaio ^ 3 Lq/Mq. We find (cf. Fig. 0) that our 
short lifetime model with tq — 10^'^ yr nearly reaches this 
limit; models with shorter lifetimes would require unreal- 
istically large L/Mhaio ratios. 

Upper limit on tq. Long lifetimes, on the other hand, 
require the ratio L/Mhaio to be small; this can lead to 
unrealistically large halo masses. The brightest quasars 
detected at redshifts z fa 2 — 3 have luminosities as large 
as L « 10"Lq (Pei 1995). We find (cf. Fig. 0) that in 
order to avoid the host halo masses of these bright quasars 
to exceed ^ 10^^ Mq, the lifetime cannot be longer than 
~ 10* yr. An alternative, standard argument goes as fol- 
lows. The black hole mass grows during the quasar phase 
as e*"/*^ where i^; = 4 x 10''(e/0.1) yr is the Eddington 
time. Assuming a conservative initial black hole mass of 
IMq, a lifetime longer than 10^ yr would give final black 
hole masses that are unacceptably large. 



3. THE CLUSTERING OF QUASARS 

As demonstrated in the previous section, equally good 
fits can be obtained to the luminosity function of quasars, 
assuming either a short or a long lifetime, and a power- 
law dependence of the mean quasar luminosity on the halo 
mass L cx M" . In this section, we derive the clustering of 
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quasars in our models, and demonstrate that they depend 
significantly on the assumed lifetime. 

The halos are a biased tracer of the underlying mass 
distribution, customarily expressed by Phaio(fc) = b^P{k) 
where Phaio and P are the halo and mass power spectra as 
a function of wavenumber k. The bias parameter for halos 
of a given mass M at a given redshift z is given by (Mo & 
White 1996) 



b{M, z)^l + 



1 



Sr. 



a{M)D{z) 



- 1 



(2) 



where D{z) is the linear growth function, a{M) is the 
r.m.s. mass fluctuation on mass-scale M (using the power 
spectrum of Eisenstein & Hu 1999), and 5c ~ 1-68 is the 
usual critical overdensity in the Press-Schechter formal- 
ism (see Jing 1999 and Sheth & Tormen 1999 for more 
accurate expressions for 6(M, z) for low M, which we find 
not to affect our results here). The bias associated with 
quasars with luminosity L in our models is given by av- 
eraging over halos of different masses associated with this 
luminosity. Following equation we obtain 



b{L,z) 



r dM^{M, z) 
Jo dM^ ' 



(3) 



dg , 



&(M,z)^(i,LM,.)/on(Af,2). 



We show in Figure |^ the resulting biasparameter b{L, z) 
in the models corresponding to Figure |l|, with short and 
long lifetimes, and at redshifts z = 2 and 3. As expected, 
quasars are more highly biased in the long lifetime model, 
by a ratio 6(long)/6(short) ^ 2. In the ACDM case, at the 
detection threshold of SDSS, we find 6(long)R:! 3 at z = 3 
and 6(long)« 2 at z = 2. Bright quasars with g' ~ 17 are 
predicted to have a bias at z = 3 as large as 6 = 10. For 
reference, we also show in this figure the bias parameters 
obtained in the OCDM cosmology, which are significantly 
lower than in the ACDM case. The number of quasars ob- 
served at a fixed flux implies an intrinsically larger number 
of sources if OCDM is assumed, because the volume per 
unit redshift and solid angle in an open universe is smaller. 
This lowers the corresponding halo mass and therefore the 
bias. 

4. COMPARISON WITH AVAILABLE DATA 

As emphasized in § ^ the presently available data leave 
considerable uncertainties in the clustering of quasars. 
Nevertheless, it is interesting to contrast the results of the 
previous section with preliminary results from the already 
relatively large, homogeneous sample of high-redshift 
quasars in the 2dF survey (Boyle et al. 2000). Our pre- 
dictions are obtained by relating the apparent magnitude 
limit to a minimum absolute luminosity at a given redshift: 
log[L„,in(z)/LB,0] = 0.4[5.48 - B + 5 log(dL(^)/pc) - 5]. 
This relation assumes no K-correction, justifled by the 
nearly flat quasar spectra {uFu — const) at the relevant 
wavelengths (e.g. Elvis et al. 1994; Pei 1995). In our 
model, the correlation length rg is given implicitly by 



e,(ro) = ¥{z)D\z)U{ro) = 



(4) 



where £,rn{r) is the usual dark matter correlation func- 
tion, and b{z) is the value of the bias parameter b[L, z) as 
determined in the previous section, but now averaged over 
all quasars with magnitudes brighter than the detection 
limit. 



b{z) 



dL- 



dL 



dL^b{L,z). (5) 



In Figure ^ we show the resulting correlation lengths 
in the long and short lifetime models. Also shown is a 
preliminary data-point with la error-bars from the 2dF 
survey, based on sa3000 quasars with apparent magnitudes 
B < 20.85 (Croom et al. 1999). The upper panel shows 
the results in our flducial ACDM model with predictions 
for this magnitude cut. The published results for rg are 
cosmology dependent, and we have simply converted them 
for our cosmological models by taking the correspond- 
ing average of the redshift-distance and angular-diameter- 
distance - this crude treatment is adequate given the large 
measurement errors. Our models generically predict a 
gradual increase of the correlation length with redshift 
("positive evolution"). The clustering is dominated by 
the faintest quasars near the threshold luminosity; as a 
result, the fixed magnitude-cut of i? = 20.85 corresponds 
to more massive, and more highly clustered halos at higher 
redshifts. There are two additional effects that determine 
the redshift-evolution of clustering: (1) quasars of a fixed 
luminosity are more abundant towards high-z, requiring 
smaller halo masses to match their number density, and 
(2) halos with a fixed mass are more highly clustered to- 
wards high-z. We find, however, that these effects are 
less important than the increase in Mhaio caused by fixing 
the apparent magnitude threshold, which gives rise to the 
overall positive evolution. 

The clustering in the long-lifetime model is stronger, 
and evolves more rapidly than in the short-lifetime case. 
As we can see, the present measurement error-bars are 
large - the whole range of life-times from 10^'^ to 10® yr 
is broadly consistent with the data, to within ^ 2a. For 
the ACDM model, the 2dF data-point is consistent with 
a fifetime of around 10^ '^±° ® yr; in the OCDM model the 
lifetime is somewhat higher, 10®^'' ® yr. It is worth empha- 
sizing here that the constraints on lifetime from clustering 
depends on the underlying cosmological parameters one 
assumes, which future large scale structure measurements 
(from e.g. the microwave background, galaxy surveys and 
the Lyman-a forest) will hopefully pin down to the accu- 
racy required here. 

Lastly, we note that observational results on tq are com- 
monly obtained by a fit to the two-point correlation of the 
form ^(r) = (r/rp)"'*'. Since ro is the correlation length 
where the correlation is unity, we expect our formalism to 
begin to fail on such a scale, because neither linear fluctu- 
ation growth nor linear biasing holds. On the other hand, 
it is also unclear whether rg as presently measured from a 
2-parameter fit to the still rather noisy observed two-point 
correlation necessarily corresponds to the true correlation 
length. While our crude comparison with existing data 
in Figure ^ suffices given the large measurement errors, 
superior data in the near future demand a more refined 
treatment, which is the subject of the next section. 
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5. EXPECTATIONS FROM THE SDSS AND 2DF 

Although existing quasar clustering measurements still 
allow a wide range of quasar lifetimes, and do not pro- 
vide tight constraints on our models, forthcoming large 
quasar samples from SDSS or the complete 2dF survey 
are ideally suited for this purpose. Here we estimate the 
statistical uncertainties on the derived lifetimes, using the 
three-dimensional quasar power spectrum Pgik) derived 
from these surveys. 

The variance of the power spectrum is computed by 



(6) 



where the large scale fluctuations are assumed to be 
Gaussian, nj- is the number of independent modes, n is 
the mean number density of observed quasars, Pqik) — 
rP{k) is the quasar power spectrum and P{k) is the mass 
power spectrum (Feldman, Kaiser & Peacock 1994). For 
a survey of volume V, and a fc-bin of size Afc, we use 
rifc = fc^Afcy/47r^. The fractional variance is therefore 
{1 + In terms of minimizing this error. 



increasing the luminosity cut of a survey has the advantage 
of raising the bias 6, but has the disadvantage of decreasing 
the abundance n. In practice, h changes relatively slowly 
with mass (slower than h cx M) whereas n varies with mass 
much more rapidly [n ~ 1/M, or steeper if M > A/*). As 
a result, we find that for our purpose of determining the 
clustering and the quasar lifetime, it is better to include 
more (fainter) quasars. 

The power spectrum Pqik) of quasars in ACDM is 
shown in Figure |^ at two different redshifts near the peak 
of the comoving quasar abundance, z — 2 and z = 3. 
We assume that redshift slices are taken centered at each 
redshift with a width of Az = 0.5 (which enters into the 
volume V above) . Results are shown in the long and short 
lifetime models, together with the expected la error bars 
from SDSS (crosses). Also shown in the lower panel are 
the expected error bars from 2dF (open squares) , which are 
slightly larger because of the smaller volume (for SDSS, we 
assume an angular coverage of tt steradians, and for 2dF, 
an area of 0.23 steradians). We do not show error bars for 
2dF beyond k ~ 0.01 Mpc/h because larger scales would 
likely be affected by the survey window. We also only 
show scales where the mass power spectrum and biasing 
are believed to be linear. 

As these figures show, the long and short tq models 
are easily distinguishable with the expected uncertainties 
both from the SDSS and the 2dF data, out to a scale of 
~ lOOMpc. The measurement errors at different scales 
are independent (under the Gaussian assumption), and 
hence, when combined, give powerful constraints e.g. for- 
mally, even models with lifetimes differing by a few percent 
can be distinguished with high confidence using the SDSS. 
However, systematic errors due to the theoretical model- 
ing are expected to be important at this level, which we 
will discuss in §^ 

In Figure |^, we have used the magnitude cuts for spec- 
troscopy, i.e. B = 20.85 for 2dF and B = 20 A [g' w 19) for 
SDSS. If photometric redshifts of quasars are sufhciently 
accurate, the magnitude cuts can be pushed fainter, fur- 
ther decreasing the error-bars - although it is likely that 

^see http:/ /asc. harvard. cdij and http://xmm.vilspa.csa.ef, respectively 



the photometric redshift errors will be large enough that 
one can only measure the clustering in projection, in some 
well-defined redshift-bin picked out using color informa- 
tion. Color selection of z > 3 quasars has already proven 
to be highly effective (Fan et al. 1999); a high redshift 
sample can give valuable information on the evolution of 
the quasar clustering (see Fig. ||). 

6. QUASAR CLUSTERING IN X-RAY 

Both the luminosity function (e.g. Miyaji et al. 2000), 
and the clustering of quasars (e.g. Carrera et al. 1998) 
has been studied in the X-ray band, analogously to the 
optical observations described above. At present, the ac- 
curacy of both quantities are inferior to that in the optical. 
Nevertheless, it is interesting to consider the clustering of 
quasars in the X-ray regime, because (1) applying the ex- 
ercise outlined above to a different wavelength band pro- 
vides a useful consistency check on our results, (2) obser- 
vations with the Chandra X-ray Observatory (CXO) and 
XMM can potentially'^ probe quasars at redshifts higher 
than currently reached in the optical, and (3) X-ray ob- 
servations are free from complications due to dust extinc- 
tion, although other forms of obscuration are possible (see 
§0). In addition, we will consider here the auto-correlation 
of the soft X-ray background (XRB) as another potential 
probe of quasar clustering and lifetime. 

The formalism we presented in § ^ and § ^ is quite gen- 
eral, and we here apply it to the soft X-ray luminosity 
function (XRLF) from Miyaji et al. (2000). The details of 
the fitting procedure are given in the Appendix. In anal- 
ogy with the optical case, we find that the clustering of 
X-ray selected quasars depends strongly on the lifetime. 
As an example, including all quasars whose observed flux 
is above 3 x 10"^'' ergcm"^ s~^, we find the correlation 
length at z = 2 to be « 4/i~^Mpc in the short lifetime, 
and « ll/i~^Mpc in the long lifetime case (ACDM). Cur- 
rent data probes the clustering of X-ray quasars only at 
low redshifts (z ^ 1, see Carrera et al. 1998), where our 
models suffer from significant uncertainties due to the sub- 
halo problem discussed in §|^. Constraints at z ^ 2 could 
be available in the future from CXO and XMM, provided 
that a large area of the sky is surveyed at the improved 
sensitivities of these instruments. 

We next focus on the quasar contribution to the soft 
X-ray background and its auto-correlation, which as we 
will see is dominated by quasar contributions at some- 
what higher redshifts. The mean comoving emissivity at 
energy E from all quasars at redshift z, typically in units 
of keVcm^^ s^^ sr~^, is given in our models by 



dL^LxiE,L), 



(7) 



where Lx(E,L) is the luminosity (in keVs~^keV~^) at 
the energy i5 of a quasar whose luminosity at (1 -I- z) keV 
is L. We have used the mean spectrum of Elvis et al. 
(1994) to include a small K-correction when computing 
the background at observed energies E ^ IkeV. The mean 
background is the integral of the emissivity over redshift. 



dx 



(1 + ^) 



(8) 
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where I is typically given in units of keV cm^^ s^^ sr^^ keV" 
Ez — E{1 + z), and x is the comoving distance along the 
line of sight. 

If S{z) is the mass fluctuation at some position at red- 
shift z, then the fluctuation of the emissivity at the same 
position and redshift is given by bx{z)S{z)j{Ez, z), where 
we have defined the X-ray emission-weighted bias bx{z) 
as 

~bx{z) = ^—-— dL-^Lx{E,,L)bx{L,z). (9) 

For simplicity, we compute the auto-correlation wg of 
the XRB using the Limber approximation, together with 
the small angle approximation, as: 

C,{E)^I{Ey' I ^W^iE^,x)Po{eK) (10) 

wg{E)= J ^-^CiiE)Joi£e), 

where is the angular power spectrum, Jq is the ze- 
roth order Bessel function, Po(^/?'x) is the linear mass 
power spectrum today, is the angular diameter dis- 
tance (= X for a flat universe), and W{E^,x) — 
j[Ez, z)bx{z)D{z) / {1 + z). When the power spectrum is 
measured in practice, shot-noise has to be subtracted or 
should be included in the theoretical prediction, whereas 
the same is not necessary for the angular correlation ex- 
cept at zero-lag. 

Our models predicts the correct mean background spec- 
trum /(-E), computed from equation ^, at i5 = 1 keV. We 
have included all quasars down to the observed 1 keV flux 
of 2 X 10~^^ erg cm~^ s~^, i.e. we used our models to 
extrapolate the XRLF to two orders of magnitude fainter 
than the ROSAT detection threshold for discrete sources 
(Hasinger & Zamorani 1997), to make up the remaining 
~ 50% of the XRB at IkeV. Our models predict a faint- 
end slope that is steeper than the Miyaji et al. (2000) 
fitting formulae, allowing faint quasars to contribute half 
of the background. The emissivities peak at z « 2, coincid- 
ing with the peak of the XRLF, implying that our model 
produces most of the XRB, as well as its auto-correlation 
signal at z « 2. Note that the known contribution from 
nearby galaxy clusters is ^ 10% (Gilli et al. 1999), which 
we ignore here. 

In Figure ^, we show our predictions for the two point 
angular correlation wg of the XRB from quasars at 1 keV. 
Most measurements at the soft X-ray bands have yielded 
only upper limits, which are consistent with our predic- 
tions (e.g. variance at ^ 0.12 at a scale of 10 arcmin. and 
E = 0.9 — 2keV, from Carrera et al. 1998; see also refer- 
ences therein). Soltan et al. (1999) obtained angular corre- 
lations significantly higher than previous results (dashed 
curve), which, taken at face value, would imply quasar 
lifetimes tq ^ 10^ yr. However, the results of Soltan et 
al. (1999) could be partially explained by galactic con- 
tamination (Barcons et al. 2000). We therefore view this 
measurement as an upper limit, which is consistent with 
models using both lifetimes we considered. Figure || shows 
that Wg predicted in the long and short lifetime models dif- 
fer by a factor of ~ 2 on angular scales of 0.1-1 degrees. 



offering another potential probe of quasar lifetimes, pro- 
vided that We can be measured more accurately in the 
future, and that the contribution to the clustering signal 
from nearby non-quasar sources (e.g. clusters) is small or 
can be subtracted out. Finally, we note that there have 
been detections of of clustering on several-degree-scales at 
the hard X-ray bands (2-10 keV) from the HE AO satellite 
(Treyer et al. 1998) - while a prediction for such energies 
would be interesting (see also Lahav et al. 1997), it would 
require an extrapolation of the X-ray spectrum, since we 
normalize by fitting to the soft-Xray luminosity function. 

7. FURTHER CONSIDERATIONS 

We have shown above that the quasar lifetime could be 
measured to high precision, using the soon available large 
samples of quasars at z ^ 3; either from the 2dF or the 
SDSS survey. This precision, however, reflects only the 
statistical errors in the simple model we have adopted for 
relating quasars to dark matter halos. The main hindrance 
in determining the quasar lifetime will likely be systematic 
errors; here we discuss how several potential complications 
could affect the derived lifetime. 

Obscured sources. Considerations of the hard X-ray 
background have led several authors to suggest the pres- 
ence of a large population of "absorbed" quasars, neces- 
sary to fit the hard slope and overall amplitude of the 
background. Although not a unique explanation for the 
XRB, this would imply that the true number of quasars 
near the faint end of both the optical and soft X-ray LF 
is ~ 10 times larger than what is observed; 90% being 
undetected due to large absorbing columns of dust in the 
optical, and neutral hydrogen in the soft X-rays (see, e.g. 
Gilli et al. 1999). Unless the optically bright and dust- 
obscured phases occur within the same object (Fabian & 
Iwasawa 1999), this increase would have a direct effect on 
our results, since we would then need to adjust our fitting 
parameters to match a ~ 10 times higher quasar abun- 
dance. We find that this is easily achieved by leaving xq 
and a unchanged, and instead raising the lifetime from 
10^'^ to 10^'^ yr in the short lifetime model, and from 10* 
to 10*'^ yr in the long lifetime model. In the latter case, 
a 10-fold increase in the duty cycle requires only an in- 
crease in tq by a factor of « 4, owing to the shape of 
the age distribution dpa/dt (Lacey & Cole 1993). Our re- 
sults would then hold as before, but they would describe 
the two cases of tq = 10^-^ and 10*'^ yr. Interestingly, 
this scenario would imply that the quasar lifetime could 
not be shorter than tq w 10^'^ yr, simply based on the 
abundance of quasars (cf. § ^). Future infrared and hard 
X-ray observations should help constrain the abundance of 
obscured sources, and reduce this systematic uncertainty. 

Multiple BH's in a single halo. A possibility that could 
modify the simple picture adopted above is that a single 
halo might host several quasar black holes. A massive (e.g. 
lO^"* Mq) halo corresponds to a cluster of galaxies; while 
the Press-Schechter formalism counts this halo as a sin- 
gle object. If quasar activity is triggered by galaxy-galaxy 
mergers, a massive Press-Schechter halo, known to contain 
several galaxies, could equally well host several quasars 
(e.g. Cavaliere & Vittorini 1998). There is some obser- 
vational evidence of perhaps merger driven double quasar 
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activity (Owen et al. 1985; Comins & Owen 1991). One 
could therefore envision that quasars reside in the sub- 
halos of massive "parent" halos - a scenario that would 
modify the predicted clustering. To address these issues in 
detail, one needs to know the mass-function of sub-halos 
within a given parent halo, as well as the rate at which 
they merge and turn on. In principle, this information 
can be extracted from Monte Carlo realizations of the for- 
mation history of halos in the extended Press-Schechter 
formalism (i.e. the so called merger tree) together with 
some estimate of the time-scale for mergers of sub-halos 
based on, for instance, dynamical friction (e.g. Kauffmann 
& Haehnelt 2000). Here, we consider two toy models that 
we hope can bracket the plausible range of clustering pre- 
dictions. 

To simplify matters, we ignore the scatter in L-M in 
the following discussion. Suppose one is interested in 
quasars of a luminosity L at redshift zq, which corre- 
spond to Press-Schechter halos of mass Mq in our for- 
malism as laid out in §|[ This choice of Mq matches 
the abundance of quasars, expressed approximately as 
L{d<^/dL) w Mo[dA^(Mo,zo)/dMo](tQ/iHub) (the merger, 
or activation rate of halos is approximated as ~ ^Hub' 
where inub is the Hubble time), and implies the bias 

6o(L)«6(Afo,^o). 

In model A, we suppose the Press-Schechter halos are 
identified at some earlier redshift zi - these would be 
sub-halos of those Press-Schechter halos identified at zq. 
Quasars of luminosity L now correspond to sub-halos of 
mass Ml . The abundance of these sub-halos is given by the 
Press-Schechter mass function dN{Mi, zi)/dMi, which is 
related to the mass function at zq by dN{Mi, zi)/dMi = 
J^^ dM [dN{M, zo ) /dM] [dJ\f (Ml , zi I M, zo ) /dMi ] where 
dMi X dA/'(Mi,zi|M, zo)/dMi is the average number of 
Ml ± dMi/2 sub-halos within parent halos of mass AI, 
given by (e.g. Sheth & Lcmson 1999) 

dA/'(Mi,zi|M,zo)/dMi = (11) 
Ml Si-6o 
Ml V2^[a2(Mi)-a2(M)]3/2 

f (^i-^o)^ ] da^{M,) 

^^"^X 2[a2(Mi) -cr2(M)] J dMi ' 

where 6i = Sc/D{zi) and So — 6c/D{zo) (see cq. ^). 

To match the abundance of quasars at luminosity L, 
we impose the condition that Mi[dN{AIi, zi) /dMi] « 
Mo[dN{MQ, zo)/dMo] , which determines Mi given zi, Mq 
and Zq. The bias of the quasars is no longer bo{L) « 
6(Mo, Zq), but is instead given by 



necessarily larger than our original bias bo{L), despite the 
fact that the bias of sub-halos should be boosted by their 
taking residence in bigger halos. This is because the rele- 
vant masses here (e.g. Mq) are generally large, and we 
find that the number of halos of mass Mq at zi is al- 
ways smaller than the number of halos of the same mass 
at Zq < zi. Hence, to match the observed abundance of 
quasars at the same L, Mi must be chosen to be smaller 
than Mq. As Fig. ^ shows, this could, in some cases, more 
than compensate the increase in clustering due to mas- 
sive parent halos. Because of these two opposing effects, 
the bias does not change by more than about 50% even 
if one considers zi as high as 10. This translates into a 
factor of ~ 2 uncertainty in our predictions for the quasar 
power spectrum. Our clustering predictions for the short 
and long lifetime models differ by about a factor of ~ 5, 
implying that the lifetime can still be usefully constrained 
at z ^ 2. As we can see from Fig. ^ at lower redshifts, or 
equivalently lower Mq/M*, our predictions for the quasar 
power spectrum should be more uncertain. 

One might imagine modifying the above model by al- 
lowing mergers to take place preferentially in massive 
parents, and therefore boosting the predicted bias. In 
model B, we adopt a more general procedure of match- 
ing the observed quasar abundance by L{d^/dL) = 
Ml /j^^ dM[dN{M, Zq) /dM] [dAf (Mi , zi |M, zq) /dMi] 
(iQ/iHub)/(Mi, M) where (iQ/tHub)/(Mi, M) is the prob- 
ability that an Mi sub-halo residing within a parent halo 
of Mq harbors an active quasar of luminosity L. It is 
conceivable that / increases with the parent mass Mq - 
a more massive parent might encourage more quasar ac- 
tivity by having a higher fraction of mergers or collisions. 
The following heuristic argument shows that one might 
expect [d7V(Mi, zi|M, zq) / dMi\f{Mi, M) to scale approx- 
imately as ~ M^/^. Let Nh be the number of sub-halos 
inside a parent halo of mass M . The rate of collisions 
is given by N^Vhcrh/R^ where Vh is the velocity of the 
sub-halos, ah is their cross-section and is the volume 
of the parent halo. Using Nh oc M (which can be obtained 
from eq. ^ in the large M limit), Vh oc a/ M/R (virial 
theorem) and (x M (fixed overdensity of ^ 200 at the 
redshift of formation) , the rate of collisions scales with the 
parent mass as M*'''^, if one ignores the possibility that 
Uh might depend on the parent mass as well. A similar 
scaling of M^ '^ has been observed in simulations of the 
star-burst model for Lyman-break objects (Kolatt et al. 
1999, Weschler et al. 1999). 

To model the enhanced rate of collisions inside mas- 
sive parent halos, we can simply modify model A by using 
/(Mi,M) = (M/Mi)i/3. The effective bias is given by 



(i,^i) 



/ dM 

J Ml 



dN{Mi,zi) 



dMi 

dN{M, Zq) dAf (Ml , zi |M, Zo) 



(12) 



dM 



dMi 



KM,zo). 



We show in Fig. ^ the ratio of bfg{L, zi)/bo{L) as a 
function of zi, for zq = 3 and zq = 2 respectively, and 
for a range of masses Mq which are representative of the 
halos that dominate our clustering signal in previous dis- 
cussions. It is interesting how the bias 6^(L, zi) is not 



b^s{L,zi) 



°° dM ^^^^^' dAA(Mi, zi|M, ZQ 
Ml dM dMi 



X / dM 

'Ml 



dM 

dN{M, zo) dAf{Mi,zi\M, zq] 



dM 



dMi 



M 



(13) 



^) KM,zq) 



We find that the above prescription does not signifi- 
cantly alter our conclusions following from model A: the 
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M^/'^ enhancement of the activation rate inside massive 
parent halos turns out to be relatively shallow, and trans- 
late to a small effect in the bias. Finally, we note that zi 
above could in principle depend on A/q and Mi , a possibil- 
ity that would require further modeling and is not pursued 
here. 

Galaxies without BH's. Another possibility that could 
modify our picture is that only a fraction / < 1 of the 
halos harbor BH's; the duty cycle could then reflect this 
fraction, rather than the lifetime of quasars. Although 
there is evidence (e.g. Magorrian et al. 1998) that most 
nearby galaxies harbor a central BH, this is not necessarily 
the case at redshifts z — 2 — 3: the fraction / of galaxies 
hosting BH's at z = 2 — 3 could, in principle have merged 
with the fraction 1 — / of galaxies without BH's, satisfying 
the local constraint. 

Using the extended Press-Schechter formalism (Lacey 
& Cole 1993), one can compute the rate of mergers be- 
tween halos of various masses. On galaxy-mass scales, the 
merger rates at z = 2 — 3 are comparable to the recipro- 
cal of the Hubble time (cf. Fig. 5 in Haiman & Menou 
2000), implying that a typical galaxy did not go through 
numerous major mergers between z = 2 — 3 and z = 0, i.e. 
that the fraction / cannot be significantly less than unity 
at z = 2 — 3. A more detailed consideration of this issue 
is beyond the scope of this paper; we simply note that the 
lifetimes derived here scale approximately as 1//, where / 
is likely of order unity. 

Larger scatter in L/M. The scatter a we assumed 
around the mean relation between quasar luminosity and 
halo mass is motivated by the scatter found empirically for 
the Mbh — Mbuigo relation (Magorrian et al. 1998). It is 
interesting to consider the sensitivity of our conclusions 
to an increased a. In general, scatter raises the num- 
ber of quasars predicted by our models, by an amount 
that depends on the slope of the underlying mass function 
dN/dM. As a result, increasing a raises, and flattens the 
predicted LF. We find that an increase of a from 0.5 to 1 
(an additional order of magnitude of scatter) can be com- 
pensated by a steeper Lm relation, typically replacing a 
with w a — 0.5. As a result of the increase in cr, quasars 
with a fixed L are, on average, associated with larger, and 
more highly biased halos. 

Nevertheless, we find that the mean bias b of all sources 
above a fixed flux (cf. eq. P), and therefore the correlation 
length To, is unchanged by the increased scatter (at the 
level of ~ 3%). The reason for the insensitivity of ro to 
the amplitude of the scatter can be understood as follows. 
The mean bias b of all sources with L > Lmin is dominated 
by the bias b{L) of sources near the threshold L,„in. The 
latter is obtained by averaging b{M) over halos of different 
masses (cf. eq. ^ , and it is dominated by the bias of the 
smallest halos within the width of the scatter, i.e. of ha- 
los with mass M^i^ ~ M/10°', where M defines thejnean 
relation between and halo mass i.e. Lmin = L{M), 

and cr quantifies the scatter (see eq. 16). As mentioned 
before, increasing the scatter makes the luminosity func- 
tion flatter, which means to match the observed abundance 
of halos at a fixed luminosity £min, one has to choose a 
higher M. In other words, M scales up with the scatter, 
and it turns out to scale up approximately as 10'^, making 



Mmin and hence the effective bias roughly independent of 
scatter. 

We note that the relation between quasar luminosity 
and halo mass can, in principle be derived from observa- 
tions, by measuring Mhaio for the hosts of quasars (e.g. by 
weak lensing, or by finding test particles around quasars, 
such as nearby satellite galaxies). 

Mass and redshift dependent lifetime. In all of the above, 
we have assumed that the quasar lifetime is a single pa- 
rameter, independent of the halo mass. This is not un- 
reasonable if the Eddington time, the timescale for the 
growth of black hole mass, is indeed the relevant time- 
scale, 4 X 10^ (e/0.1) yr. Implicit in such reasoning is that 
the active phase of the quasar is coincident with the phase 
where the black hole gains most of its mass. This is not 
the only possibility; sec Hachnelt et al. (1999) for more 
discussions. One can attempt to explore how tg depends 
on halo mass by applying our method to quasars grouped 
into different absolute luminosity ranges, but the intrinsic 
scatter in the mass-luminosity relation should be kept in 
mind. We emphasize, however, since we fit the luminosity 
function and clustering data at the same redshift, there 
is no need within our formalism to assume a redshift in- 
dependent lifetime. In fact, performing our exercise as a 
function of redshift could give interesting constraints on 
how tg evolves with redshift. 

8. CONCLUSIONS 

In this paper, we have modeled the quasar luminosity 
function in detail in the optical and X-ray bands using 
the Press-Schechter formalism. The lifetime of quasars tg 
enters into our analysis through the duty-cycle of quasars, 
and we find that matching the observed quasar LF to dark 
matter halos yields the constraint 10^ ^ tg <, 10^ yr: 
smaller lifetimes would imply overly massive BH's, while 
longer lifetimes would necessitate overly massive halos. 
This range reassuringly brackets the Eddington timescale 
of 4 X 10^ (e/0.1) yr. 

The main conclusion of this paper is that the life- 
time (and hence e, if the Eddington time is the relevant 
timescale for quasar activity) can be further constrained 
within this range using the clustering of quasars: for 
quasars with a fixed luminosity, longer tg implies larger 
host halo masses, and higher bias. We find that as a result, 
the correlation length tq varies strongly with the assumed 
lifetime. Preliminary data from the 2dF survey already 
sets mild constraints on the lifetime. Depending on the 
assumed cosmology, we find tg = IQ'^-'^^O-^ yr (ACDM) 



or tc 



yr (OCDM) to within la statistical un- 
certainty. These values are also found to satisfy upper 
limits on the auto-correlation function of the soft X-ray 
background. 

Forthcoming large quasar samples from SDSS or the 
complete 2dF survey are ideally suited for a study of 
quasar clustering, and they can, in principle constrain the 
quasar lifetime to high accuracy, with small statistical er- 
rors. We expect the modeling of the quasar-halo relation, 
as well as the possible presence of obscured quasars, to 
be the dominant sources of systematic uncertainty. Not 
discussed in depth in this paper is the possibility of using 
higher moments (such as skewness), which our models also 
make definite predictions for, and will be considered in a 
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future publication. Remarkably, our best determination of 
the lifetime of quasars might come from the statistics of 
high-redshift quasars, rather than the study of individual 
objects. 

Near the completion of this work, we became aware of 
a similar, independent study by P. Martini & D. Wein- 
berg. We thank M. Haehnelt, K. Menou, E. Quataert, 
U-L. Pen., U. Seljak, R. Sheth and D. Spergel for useful 
discussions, T. Miyaji for his advice on the X-ray luminos- 



ity function, and T. Shanks and S. Groom for discussions 
on the 2dF survey. Support for this work was provided by 
the DOE and the NASA grant NAG 5-7092 at Fermilab, 
by the NSF grant PHY-9513835, by the Taphn Fellowship 
to LH and by NASA through the Hubble Fellowship grant 
HF-01119.01-99A to ZH, awarded by the Space Telescope 
Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc., for NASA 
under contract NAS 5-26555. 
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APPENDIX 

In this Appendix, we describe our models for the luminosity function (LF) of quasars, based on associating quasar 
BHs with dark halos. Our main assumption is that there is, on average, a direct monotonic relation between halo 
mass and quasar light. Our treatment is similar to previous works (Haiman & Loeb 1998; Haehnelt et al. 1998), but 
differs in some of the details. We adopt the parameterization of the observational LF in the optical B band, given in 
the redshift range < z ^ 4 by Pei (1995). We assume the background cosmology to be either flat (AGDM) with 
(rJA,nni,/i,cr8h-i,n) = (0.7,0.3,0.65,1.0,1.0) or open (OGDM) with (f^A, f^m, /i, crgh-i , n) = (0,0.3,0.65,0.82,1.3). The 
LF quoted by Pei (1995) is scaled appropriately with cosmology by keeping {d(t)/dL)dV(if^=const, where dV is the volume 
element, and c^l is the luminosity distance), so that {d(j)/dL)dL is the comoving abundance in Mpc~^ of quasars with 
B-band luminosity L (in solar units Lb,q)- The comoving abundance dN/dM{M, z) of dark halos is assumed to follow 
the Press-Schechter (1974) formalism. We assume that each halo harbors a single quasar that turns on when the halo 
forms, i.e. triggered by merger (e.g. Percival & Miller 1999), and shines for a fixed lifetime tg (relaxing these assumptions 
is discussed below in section The duty-cycle /on of halos with mass M at redshift z, is then given by the fraction of 
these halos younger than tg. The distribution of ages dpa/dt{M, z, t) for halos of mass M existing at redshift z is obtained 
using the extended Press-Schcchter formalism, which assumes that the halo formed at the epoch when it acquired half of 
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its present mass (Lacey & Cole 1993). The duty-cycle, which is the probability that a dark matter halo of a given mass 
harbors an active quasar, is simply 

/on(M,z)- dt^{M,Z,t). (14) 

A model in which the quasar turns on/ofF more gradually (as expected if the mass of the BH grows significantly during 
the luminous quasar phase) is equivalent to one having additional scatter in the ratio £/Af, which is discussed in § |^. We 
next relate the quasar luminosity to the mass of its host halo. We define dp/dL{M, L, z) to be the probability that a halo 
of mass M at redshift z hosts a quasar with luminosity L, and express this quantity as 

^{L,M,z) = ^(L,Lm,.)/o„(M,z). (15) 

Here dg/dL{L, Lm.z) is the probability distribution of luminosities associated with the subset of halos of mass M 
harboring a live quasar (normalized to dLdg/dL — 1), and Lm,z is the mean quasar luminosity for these halos. In the 
limit of a perfect intrinsic correlation, we would have dg / dL{L, Lm,z) = ^{L — Lm,z)', more realistically, this correlation 
will have non-negligible scatter. Lacking an a-priori theory for this scatter, we here simply assume that it follows the 
same functional form as the scatter found empirically for the Afbh — Afbuigo relation (Magorrian et al. 1998), and we set 

%{L,Lm) (X exp(-(logL/LM,.)V2fT2)- (16) 
dL 

For reference, we note that the empirical scatter between Mbh and Afbuigc gives a ^ 0.5, it is not yet clear, however, 
what fraction of this scatter is intrinsic vs. instrumental (van der Marel 1999). One might expect the scatter in the 
L-Mhaio relation not to be significantly larger, since (i) for a sufficiently high fueling rate, the luminosity L corresponding 
to Mbh is likely to always be near the Eddington limit, and (ii) at least for disk galaxies, the bulge luminosity correlates 
well with the total luminosity itot (c ~ 0.5, see e.g. Andredakis & Sanders 1994); Ltot is tightly correlated with the 
velocity dispersion through the TuUy-Fisher relation (e.g. Raychaudhury et al. 1997) as is the total halo mass to 
(Eisenstein & Loeb 1996). Nevertheless, in § |^ below we will investigate the consequences of an increased scatter. We 
note that an extension of the models presented here, by following the merger histories of halos and their BH's can, in 
principle, be used to estimate the scatter in L/Mhaio- Cattaneao, Haehnelt & Rees (1998) have used this approach to fit 
the observed relation Mbh/Mbuigo, including its scatter. 

Under the above assumptions, the cumulative probability that a halo of mass M hosts a quasar with luminosity equal 
to or greater than L is given by 

piL,M,z)^ fon{M,z) f dL^{L,LM.z), (17) 



L 



and matching the observed cumulative quasar LF requires 



L 



d(b poc c/ZV 

dL^{L,z)^J^ dM — {M,z)p{L,M,z), (18) 



or alternatively, matching the differential LF gives 

^{L, z) = dM^iM, z)^{L, Lm,.)/o„(M, z). (19) 

Equation |l^ or |l^, together with equations |lj, |l^, and |l^ implicitly determines the function Lm,z: once the quasar 
lifetime tq aSd magnitude of scatter a are specmedT^ In general, these equations would need to be solved iteratively. In 
practice, we have found that sufficiently accurate solutions (given the error bars on the observational LF in Pei 1995; cf. 
Fig. |l|) can be found by using the simple power-law ansatz: 

7i,f \ a(z) 
JWhalo 



Lm,z = xo(z)Mhaio y-^J ' (20) 

where the coefficients xo{z) and a(z) depend on tg, a, and the underlying cosmology (Haiman & Loeb 1998; Haehnelt 
et al. 1998). In summary, assuming a fixed scatter cr, our model has only one free parameter, the quasar lifetime tq. 

We emphasize that our parameterization in equation is purely phenomenological - it gives us a convenient way 
to relate the quasar luminosity to the host halo mass {Lm,z)- In reality, the quasar luminosity likely depends on the 
details of its immediate physical environment (e.g. gas supply, magnetic fields, angular momentum distribution, etc.), in 
addition to the halo mass. Our description includes these possibilities only in allowing a non-negligible scatter around the 
mean relation Lm,z- The rationale behind this choice is that the average properties of the physical environment should 
ultimately be governed by the halo mass (or circular velocity), as expected within the picture of structure formation via 
hierarchical clustering. 
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A useful check on the physical implications of equation is obtained by assuming that the luminosity Lm,z is produced 
by a BH of mass Mbh, shining at the Eddington limit LecIcI — (47rG/iTOpc/(TT)Mbh- In the mean spectrum of a sample of 
quasars with detections from radio to X-ray bands (Elvis et al. 1994), « 7% of the bolometric luminosity is emitted in 
the rest-frame B band, resulting in L — O.OTLEdd = 5 x lO'^LB_0(Mbh/MQ). Equation pO| then translates into a relation 
between the mass of a BH and its host halo, 

Mbh = 10-^-^X0 (z)Afhaio ( j • (21) 

As an example, Haehnelt et al. (1999) argue that the central BH mass is determined by a radiative feedback from 
the central BH that would unbind the disk in a dynamical time. Their derived scaling corresponds to a = 2/3 and 
xo c>c (1 + z)^^^, not far from what we find for the long-lifetime case (cf. Figure and discussion below). 

In Figure |^, we show the values of the parameters xo{z) and a{z) obtained in our models when two different quasar 
lifetimes are assumed, tq — 10^-^ (solid curves) and tq — 10* yr (dotted curves). The filled dots show the parameters 
in ACDM, and the empty dots in the OCDM cosmology. We have set the arbitrary constant Mq = 10^^ Mq in both 
cases. Note that tg determines both a and xq, and therefore the values of a and xq are correlated. In general, the 
fitting parameters show little evolution in the range 2 < z < 4, around the peak of the quasar LF. According to equation 



21], the corresponding BH masses in, e.g. a lO^^M© halo at 2 < 2; < 4 are Mbh « 4 x lO'^'Mhaio = 4 x lO^M© and 
Mbh « 2 X 10~^Mhaio = 2 X 1O''M0 in the short and long lifetime models, respectively. 

The fitting procedure described above can be repeated in the X-ray bands. We therefore fit the XRLF using equation |20| 
analogously to the optical case, except Lm,z now denotes the X-ray luminosity at 1 keV, quoted in units of erg s~^. Note 
that the XRLF in Miyaji et al. (2000) is quoted a function of luminosity at observed 1 keV, i.e. no K-correction is 
applied (alternatively, the XRLF can be interpreted as the rest-frame luminosity function of sources with an average 



intrinsic photon index of 2 
to those shown in Figure 
from that in the optical 



Figure || show the resulting fitting parameters xq and a in the ACDM cosmology, analogous 
for the optical case. It is apparent that both parameters have a somewhat behavior different 
his reflects the fact that the mean quasar spectrum must evolve with redshift, or at least is 
black- hole/halo mass dependent: if every quasar had the same spectrum, or at least a similar X-ray/optical flux ratio, 
the fitting parameters derived from the optical and X-ray LF would differ only by a constant in xq. For our purpose of 
deriving clustering, it is sufhcient to treat xq and a as phenomenological fitting parameters, and we do not address the 
physical reason for the apparent spectral evolution (see Haiman & Menou for a brief discussion). 

It is important to note that the simple power-law ansatz in equation pO| with the parameters shown in Figure ^ adequately 
fits only the faint end of the XRLF. In the optical case, the entire range of observed luminosities is well matched by our 
models (cf. Fig ||). In comparison, the well-fitted range in X-rays typically extends from the detection threshold to up 
to 2-3 orders of magnitude in luminosity (i.e. typically upto '--^ 3 x 10^^ erg s~^), depending on redshift, and our models 
underestimate the abundance of still brighter quasars. One might then consider searching for a different ansatz to replace 
equation ^ that fits the entire range of the observed LF. However, we have verified that the rare quasars with these high 
luminosities would contribute negligibly both to the clustering signals, or the XRB investigated here. Therefore, we did 
not consider further improvements over equation |2C|, since this would not change our results. 
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Fig. 1. — Fits to the quasar luminosity function at icdshifts z = 1 and 3 in our models, with two different quasar lifetimes Iq = lO^-S (solid 
curves) and Iq = 10* yr (dotted curves). Also shown arc the data, and fitting function (dashed curves) for the LF from Pei (1995). The 
quality of our fits at different redshifts or in the OCDM model are similar. The upper labels show the corresponding apparent magnitudes in 
the SDSS g' band, assuming that the intrinsic quasar spectrum is the same as the mean spectrum in the Elvis et al. (1994) quasar sample. 
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Fig. 2. — Bias parameter b{L, z) of quasars as a function of their intrinsic B-band luminosity (lower labels) or apparent SDSS g' magnitude 
(upper pLabels), in the models with two different quasar lifetimes tg = 10®'^ (solid curves) and tq = 10 yr (dotted curves) as shown in 
Figure hi. For comparison, the bias parameter is also shown in the OCDM model. Quasars are more highly biased in the long lifetime models, 
and in ACDM. 
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Fig. 3. — Correlation length ro in our models with two different quascir lifetimes tq = 10^-^ (soUd curves) and tg = 10* yr (dotted curves). 
An apparent magnitude cut of B < 20.85 was used, corresponding to the limits of the 2dF survey (Croom et al. 1999). The open square 
shows a preliminary result from 2dF. The upper panel shows our results in the ACDM model, and the lower panel in OCDM. 
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Fig. 4. — Three-dimensional power spectrum PQ{k) of quasars in ACDM at two different redshifts near the peak of the comoving quasar 
abundance. Results are shown for quasars with B < 20.4 (or g' < 19), in the long (dotted curves) and short lifetime (solid curves) models, 
together with the expected la error bars from SDSS (crosses) with an assumed area of tt steradians. The slightly lower curves in the lower 
panel refer to 2dF (open squares), with a magnitude cut of B = 20.5, and show the expected error bars from an assumed area of 0.23 
steradians. 
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Fig. 5. — Angular correlation function of the total XRB, w$ at E = 1 keV. The dashed lines in the upper right corner indicate w$ at 
E = 1.15 keV with quoted ilu uncertainties as measured from the ROSAT All Sky Survey (Soltan et al. 1999), which can be considered an 
upper limit. 
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Fig. 6. — The ratio of the effective bias 6^(L) (in a model allowing multiple quasars per halo) to the original bias bo{L) (for one quasar per 
halo) as a function of zi at fixed luminosity L, where zi is the redshift when sub-halos are identified. The squares (open for Mq = lO^-^'^ Mq 
and solid for Mq = IO^^-^Mq) show this relation for quasars at zq = 3, whereas the triangles are for zo = 2 (open for Mq = 10^-^ Mq and 
solid for Mo = IO^^Mq). The model is ACDM. See text for discussions. 
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Fig. 7.=^ Fitting parameters a{z) and xo{z) for the mean relation between B-band quasar luminosity Lb and halo mass M^aio given in 
equation M, for two different quasar lifetimes tg = 10^'^ (solid curves) and tg = 10* yr (dotted curves). J^he filled dots correspond to a 
ACDM ana the open dots to an OCDM cosmology. In all cases, we assumed a scatter with a = 0.5 (cf. eq. hq) around the mean L — A/halo 
relation. 
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Fig. 8. — Same as Figure M, except the X-ray quasar luminosity is used at IkeV. Only the ACDM case is shown. 



